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Dynamics and stability of icosahedral Fe-Pt nanoparticles 

Pawet T. Jochym,* a Jan Lazewski/* Matgorzata Sternik a and Przemyslaw Piekarz a 


The structure, dynamics and stability of Fe-Pt nanoparticles have been investigated using DFT-based techniques: total energy 
calculations and DFT molecular dynamics. The investigated systems included multi-shell and disordered nanoparticles of iron 
and platinum. The study is concerned with icosahedral particles with magic number of atoms (55): iron-terminated Fe 43 Pt 12 , 
platinum-terminated Fe 12 Pt 43 , and disordered Fe 27 Pt 28 . Additionally, the Fe 6 Pt 7 cluster has been investigated to probe behaviour 
of extremely small Fe-Pt particles. Molecular dynamics simulations have been performed for a few temperatures between T = 
150 — 1000 K. The calculations revealed high structural instability of the Fe-terminated nanoparticles and a strong stabilising 
effect of the Pt-termination in the shell-type icosahedral particles. The platinum termination prevented disordering of the particle 
even at T = 1000 K indicating very high melting temperatures of these Fe-Pt icosahedral structures. The analysis of evolution of 
the radial distribution function has shown significant tendency of Pt atoms to move to the outside layer of the particles - even in 
the platinum deficient cases. 


1 Introduction 

Nanoalloys are bi- or multi-component metallic nanoparticles 
(NPs), with the complex structures and properties, which de¬ 
pend on their size, structure, and composition. 15 The size of 
nanoalloys (1-100 nm) places them between small molecules 
and bulk crystals. They can be considered as the intermediate 
phase sharing common properties of atoms and solid materi¬ 
als. By adjusting the size, geometry and chemical composi¬ 
tion, nanoalloys can be optimised for various applications in 
catalysis,® nanomedicine, 7 optics,® data recording 9 and en¬ 
ergy storage. E® 

Iron-platinum (Fe-Pt) nanoalloys consist of two metallic 
components with significant ratios (Pt/Fe) of masses (195.1 
/ 55.8 in amu), bulk moduli (230 / 170 in GPa) and atomic 
radii (1.46 / 1.30 in A) which result in very interesting, not 
only geometrical, features. They belong to the most studied 
systems due to an extremely high uniaxial magnetocrystalline 
anisotropy (K u = 7 x 10 6 J/m 3 ), which makes Fe-Pt a promis¬ 
ing material for applications in ultra-dense magnetic record¬ 
ing media. 9 As evidenced by many studies, Fe-Pt systems ex¬ 
hibit enhanced catalytical properties, comparing to pure plat¬ 
inum surfaces. Fe-Pt nanoalloys may be used as a power¬ 
ful catalyst for the electro-oxidation of formic acid in fuel 
cells 11 13 and as highly efficient oxidation catalysts for degra¬ 
dation of organic pollutants. 14 Fe-Pt NPs are also promising 
candidates for biomedical in vivo applications (e.g. for cel¬ 
lular magnetic resonance imaging 15 ) because of their high 
stability in the presence of oxygen. 16 By coating with inor¬ 
ganic molecules, the core-shell Fe-Pt nanoparticles become 
optically detectable 17 and can be functionalised for the hyper¬ 
thermia treatment of cancer and advanced radiopharmaceuti- 
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cal applications. 18 

Due to small sizes, standard diffraction techniques cannot 
be applied for NPs and most information about their struc¬ 
tural properties comes from high-resolution transmission elec¬ 
tron microscopy.® 21 Relevant information about structural 
stability and electronic properties of the Fe-Pt nanoclusters 
have been obtained by the theoretical studies based on ei¬ 
ther semiempirical potentials or ab initio techniques. 22 27 The 
Monte Carlo simulations showed that a continuous transfor¬ 
mation between the ordered structure Llo it PA) and a dis¬ 
ordered phase in Fe-Pt NPs occures at tem peratu res lower 
than the bulk melting temperature (1572 K). 23 24 Disorder¬ 
ing processes in nanoalloys are enhanced due to finite-size 
and surface effects, e.g. strong surface segregation tendency 
of Pt atoms.®®® The structural order and magnetic properties 
of Fe-Pt nanoalloys with various morphologies and chemical 
composition have been studied previously within the large- 
scale Density Functional Theory (DFT) calculations.® 7 In the 
present study, we focus on icosahedral Fe-Pt NPs, which 
are energetically favoured over the Llo phase for small par¬ 
ticles. 192021 Application of DFT-based methods, including 
DFT molecular dynamics, provides additional, independent 
information on the behaviour of Fe-Pt NPs, which supplement 
the published studies based mainly on semiempirical poten¬ 
tials and supplies a valuable consistency check on the already 
established models of these systems. 

Our study of Fe-Pt NPs has been divided into two main 
parts. The first part is a set of total energy calculations and 
structural optimisations, designed to establish static stability 
properties of the investigated systems. The methods employed 
in this part are described in section [2T] and the results in sec¬ 
tion HU The second part of the work concerns the investi¬ 
gation of dynamical stability of the systems using DFT-based 
molecular dynamics (MD). Description of the methods em- 







ployed for this part of the work is included in section [Z2| and 
the results are discussed in section [T2l Section [4] contains the 
results of the investigation of dynamical stability of the disor¬ 
dered Fe 6 Pt 7 and Fe 27 Pt 28 particles. The concluding section]?] 
summarizes the key results of the paper. 

2 Methods 

A common denominator of the methods used in the present 
work is the Density Functional Theory. 28 29 Both parts of our 
analysis use the same DFT code as implemented in Vienna Ab 
initio Simulation Package (VASP) 3031 and the same atomic 
data sets provided with this package. For the preparation and 
analysis of the molecular dynamics runs and control of the 
DFT computations, we have used standard, python-based soft¬ 
ware stack 32 ®. The details of the computational setups used 
for the work are described below - separately for each part of 
the analysis. 

2.1 Structural optimisation 

The static, zero Kelvin temperature calculations were per¬ 
formed using the full-potential projector-augmented wave 
method®® within the GGA approach in PAW-PBE form.®® 
The following valence base configurations were included: Fe 
3d 1 As 1 and Pt 6s l 5d 9 . The integrations in the reciprocal space 
was reduced to the single k-point T and the energy cut-off 
for the plane waves expansion was equal to 320 eV. The crys¬ 
tal structure was optimised using the conjugate gradient tech¬ 
nique with the energy convergence criteria set at 10 -7 and 
10 -5 eV for the electronic and ionic iterations, respectively. In 
order to prevent interactions between system images created 
due to the periodic boundary conditions used in VASP, each 
nanoparticle was placed in the big 20Ax20Ax20A box con¬ 
taining about 10 A wide vacuum. This limited force constants 
between neighbouring images below 0.3% of the strongest 
force constant in the system. 

2.2 Molecular dynamics 

The molecular dynamics simulations were carried out in the 
same DFT framework as the static calculations described 
above and with the set of methods and tools used in our previ¬ 
ous work on minerals. 41 We used the PAW-PBE atomic data 
sets and the spin-polarised VASP calculation for the DFT part. 
To lower the computational cost we adjusted the accuracy pa¬ 
rameters of the calculation - since the MD does not need ex¬ 
tremely precise inter-atomic forces. The precision parameter 
of the calculation was set to Normal, the cut-off energy was re¬ 
duced to 275 eV and the energetic convergence threshold was 
set at 10 -6 eV. A small numerical noise added by lowering 
accuracy of the calculation is drowned by the thermal noise 


of the system. Furthermore, additional energy pumped into 
the system by this noise is cancelled out by the action of the 
thermostat. The actual molecular dynamics used the VASP 
implementation of the Nose thermostat MD algorithm with 
10 fs time step, which was determined by convergence test¬ 
ing. The use of this fairly large time step was possible thanks 
to large masses of atoms in the system. The starting configu¬ 
rations correspond to perfect icosahedral geometries with the 
multi-shell or disordered structures. In case of the disordered 
structure the atomic sites were assigned randomly to one of 
the species by drawing from the pre-determined set of atoms 
- 6 Fe and 7 Pt or 27 Fe and 28 Pt, respectively. The result¬ 
ing structures were pre-optimised, by standard static methods 
(see sec. |2.1| ) to avoid large inter-atomic forces at the start of 
the MD run. For the same reasons we have used the results 
of the static structural optimisation (sec. 3.1 ) for the starting 
configuration of ordered particles. 




io 2 

10 1 

10 ° 

10' 1 

10' 2 

IO ' 3 

10' 4 


Fig. 1 Distribution of particle kinetic energy in the icosahedral-shell 
system (dots) and Maxwell-Boltzmann distribution fits for each data 
set (lines). The best-fit temperature for each data set is included in 
the legend. 

The calculation procedure involved a thermalisation run of 
the system with kinetic energies assigned to atoms accord¬ 
ing to the Maxwell distribution. To assess the level of ther¬ 
malisation of the system we calculated a distribution of ki¬ 
netic energies in each step for all the atoms and for an en¬ 
semble configurations in the 0.1 ps sliding window (10 steps). 
The average energies in each step were compared with the 
value expected for the perfect Maxwell distribution. When 
the systematic drift of the computed average value transitioned 
into random-looking fluctuations, we performed an additional 
check for good thermalisation of the system by computing the 
full distribution of kinetic energies in the ensemble built from 
the configurations in the 1 ps window (100 steps) and com¬ 
paring it with the Maxwell distribution. The comparison was 
performed by fitting the Maxwell distribution to the histogram 
of the kinetic energy in the ensemble and comparing the ob¬ 
tained temperature with the expected value. We also checked 
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Fe and Pt shells exhibit very high structural stability .^ 1 The 
stability of Ih NPs results from the low energy of the close- 
packed (111) surfaces, which compensates the internal stress 
existing in the particle core. The present study focuses on 
two types of icosahedral particles: particles with perfect shell 
structures (terminated with Pt or Fe atoms) and disordered 
NPs with approximately equal number of Fe and Pt atoms. 
The ordered Fe 12 Pt 43 particle with the alternating shell struc¬ 
ture: 1 Pt, 12 Fe and 42 Pt atoms is displayed in Fig.|4ja). 

3.1 Static calculations 


Fig. 2 Distribution of particle kinetic energy in the random systems 
(dots) and Maxwell-Boltzmann distribution fits for each data set 
(lines). ’Small’ system (left column) is a disordered Fe 6 Pt 7 cluster, 
while ’Large’ system (right column) is a disordered Fe 27 Pt 28 cluster 
(see sec. [4). The best-fit temperature for each data set is included in 
the legend. 

if the size of the residuals is consistent with the size of the 
fluctuations expected for the ensemble. 

The combination of all such comparisons is presented in 
Figs [I] and [2] The probability distributions presented in these 
figures are based on the whole measurement part of the MD 
run - the thermalisation part was discarded. Very good agree¬ 
ment of data points with the perfect Maxwell distribution is 
slightly broken only at the extreme large energy range in the 
case of the Fe-terminated system at the lowest temperature 
(T = 150 K). Notice the logarithmic scale on both axis in Figs 
[T] and [2] which was chosen to amplify any discrepancy from 
the expected distribution - on the linear scale the discrepancy 
is invisible. This is the range of a high-energy tail of the dis¬ 
tribution and one can fully expect large relative fluctuations 
there, due to the small probability density in this region and 
thus, small counts in the histogram bins. Furthermore, as the 
calculations revealed, the Fe 43 Pt 12 is a quite unstable system 
which results in difficulties with bringing it to equilibrium at 
low temperatures. Nevertheless, we do not consider this sin¬ 
gle discrepancy to be significant for the analysis carried out 
in this work and regard all investigated systems as adequately 
thermalised. 

3 Icosahedral Fe-Pt nanoparticles 

Icosahedron (Ih) is one of the five Platonic solids and one 
of the^non-crystallographic structures with fivefold symme¬ 
tries.! 42 It consists of 20 surfaces, 12 vertices and 30 edges, 
each surface has close-packed structure with the (111) ori¬ 
entation. The most stable particles are built of closed shells 

with the number of atoms N s h e ii = 1,12,42,_ Fe-Pt NPs 

crystallising in the icosahedral geometry with the alternating 


The first step, after constructing individual NPs, was to op¬ 
timise the structure at T = 0 K, by minimising the total en¬ 
ergy of the system, keeping 7} z point group symmetry ele¬ 
ments only, since 4 symmetry was broken in periodic bound¬ 
ary conditions Q The lowering of particle symmetry allows 
not only for volume changing but also for surface shell relax¬ 
ation. In both cases, Fe 43 Pt 12 and Fe 12 Pt 43 , high symmetry 
configuration with well-defined alternating Fe/Pt shells was 
obtained. The calculated radial distribution functions of opti¬ 
mised structures are presented in Fig. [3] In icosahedra there is 
only one non-equivalent position of atoms in the first shell, 
which defines the distance between the centre and the first 
shell, and two different positions in the second shell: at ver¬ 
tices and edges of the icosahedron. Interestingly, in both NPs, 
the distances from the central atom to the first shell (2.643 A 
and 2.587 A for Fe 43 Pt 12 and Fe 12 Pt 43 , respectively) as well 
as to vertices in the second shell (5.091 A for Fe 43 Pt 12 and 
5.054 A in Fe 12 Pt 43 case) are quite similar. A larger differ¬ 
ence (^ 0.1 A) is found for the distance between the centre 
and edge atoms in the second shell (4.374 A in Fe 43 Pt 12 and 
4.477 A in Fe 12 Pt 43 case). 



Fig. 3 Radial distribution of atomic positions from the mass centre 
of Fe 43 Pt 12 (left panel) andFe 12 Pt 43 (right panel). 

In the surface shell, there are two characteristic nearest 
neighbours (NN) distances: from a vertical atom to the near- 


*4 point group includes 5-fold axes, which are not compatible with the trans¬ 
lational symmetry and 3D periodic boundary conditions imposed on supercell 
in VASP package. Remaining elements of 4 group, forming T h point group, 
were considered in our static calculations. 




































est edge atoms (1-2 in Fig. [4]) and between two nearest edge 
atoms (2-2 in Fig. [ 4 ]). In Fe 43 Pt 12 , they read 2.704 A and 
2.677 A, while in Fe 12 Pt 43 - 2.767 A and 2.663 A, respec¬ 
tively. In the Pt subsurface shell NN distances are equal to 
2.778 A, while in the Fe one they are 2.720 A. For compari¬ 
son, in the bcc iron structure simulated with the same atomic 
data set, 43 the NN distance is equal to 2.485 A and in fee plat¬ 
inum 2.772 A. This means that, contrary to platinum shells, 
the inter-atomic distances in both iron shells, the surface one 
of Fe 43 Pt 12 and subsurface of Fe 12 Pt 43 , are strongly stretched 
(~ 8 — 10%) which can be one of the main reasons of the iron- 
terminated NP instability. The similar influence of strain on 
iron surface stability was observed previously in the Fe mono- 
layer on W(110) surfaced and in FeAu and FePt multilay- 
ers 45 46| j n | 30t j 1 p ar ti c i es? the nn distances between the Fe 
and Pt atoms located in the first and second shell (^ 2.50 A) 
and between the central atom and the first shell (~ 2. 60 A) 
are reduced comparing to the Fe-Pt distance in the bulk fet 
Llo structure (2.70 A).® It causes the internal stress in the 
core, which may be the driving force for the amorphisation 
processes and may induce the structural transformation ob¬ 
served in larger Fe-Pt particles. 

3.2 Molecular dynamics 

The dynamical stability of a nanoparticle is quite different 
from the static one. The later requires only that the system 
is in the local energy minimum, the former necessitates stabil¬ 
ity of the structure against finite, even large, thermal fluctua¬ 
tions. The amplitude of fluctuations is determined by the tem¬ 
perature of the system and the strength of inter-atomic bonds 
present in the system. To go beyond the static stability analysis 
described above, we need to include temperature in our calcu¬ 
lations. We used the DFT-MD technique described above (see 
section |Z2| ). The simple way to asses the effects of such cal¬ 
culation is by inspection of the evolution of the structure. 




Fig. 4 Fe 12 Pt 43 (a) and Fe 43 Pt 12 (b) structures at T = 150 K after 10 
ps thermalisation and 20 ps of MD run. The difference in structural 
stability is visible even at such a low temperature. Elements: Fe - 
dark/brown, Pt - light/grey. 1 and 2 denote two non-equivalent 
positions of atoms in the surface shell of ideal icosahedron. 

The structures after thermalisation (10 ps) and equilibration 
(20 ps) runs at target temperature T = 150 K are depicted in 
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Fig. 5 Radial distribution of atomic positions from the mass centre 
(c.m.) of the icosahedral-shell system. The results for all 
temperatures are compared side-by-side for two investigated 
systems Pt-terminated (left) and Fe-terminated (right). The 
probabilities are normalised to unity per each element. The 
semitransparent colour of the Fe plot (orange) is used to show the Pt 
plot (blue) in the same range. 



Fig-0 At the first glance the qualitative difference between 
platinum-terminated (a) and iron-terminated (b) particles is 
obvious. All the Fe atoms are in the subsurface in the for¬ 
mer, nearly perfect structure, while there are many platinum 
atoms on the surface of the later, strongly distorted structure 
- despite the fact that this is a Pt-deficient system. What is 
more, the picture is virtually identical at higher temperatures. 
Even at T = 1000 K there is no qualitative visual difference 
between Fig. [4] and the corresponding picture at that temper¬ 
ature. This calls for more detailed, quantitative analysis. In 
the crystalline material the radial distribution function (RDF) 
would be a first-line tool for such investigation. Here, the sys¬ 
tem is finite and we cannot use standard RDF directly. Instead 
we have used radial probability distribution pt(r) of finding 
atom of element i at the distance r from the centre of mass of 
the system. The probability is calculated over the ensemble of 
last 10 ps (1000 steps) of configurations at the end of the MD 
run. We performed the same calculation on the much smaller 
ensembles (0.1 — 1 ps) and determined that there is no signif¬ 
icant difference except for much larger fluctuations present in 
the histograms (for obvious statistical reasons). 

The probability density Pi(r) plots, normalized to unity for 











































each atomic species, for all investigated icosahedral systems 
and temperatures are collected in Fig. [5] It should be remem¬ 
bered that the systems are started in the perfectly ordered state 
with temperature-determined kinetic energy distribution. At 
the start of the simulation pi(r) function for all systems re¬ 
sembles the plot for Fe 12 Pt 43 at T = 150 K (top, left corner of 
Fig. [5]) or the static RDF from Fig. [3] The probability distri¬ 
butions in the Fig. [5] are calculated at the end of the equilibra¬ 
tion. Since the relaxation processes may have very different 
time scales, one cannot be sure that the true equilibrium state 
has been reached. Nevertheless, the lack of systematic drift in 
the probability distributions at the end of the simulation period 
allows us to conclude that the systems are reasonably close to 
equilibrium. 

There is a clear difference between final states of the Pt- 
terminated (left column in Fig. [5]) and Fe-terminated (right 
column in Fig. [5j system. The former, essentially, exhibits 
temperature-caused broadening of the distribution peaks, 
while the later shows clear signs of substantial reconstruction 
of the structure as well as much stronger broadening and dis¬ 
tortion of the peaks. 

Since there is not much happening in the Pt-terminated 
structure, let us concentrate on the Fe-terminated one. We 
can distinguish two main phenomena in these plots. The first 
one is a significant disordering of the system, shown in very 
broad and irregular peaks - even at low temperatures. The 
second phenomenon is a systematic drift of the platinum dis¬ 
tribution to the external shell of the particle. We expect the 
former process to have fairly short time scale and low activa¬ 
tion barriers, opposite to the other one which requires global, 
collective movements and thus large activation energy, and ex¬ 
hibits much longer time scale. Consequently, the reconstruc¬ 
tion is probably close to finished only in the high temperature 
case (T = 1000 K bottom of Fig. [5]), is under way for inter¬ 
mediate temperatures (T = 500 K), and is only starting for 
lower temperatures (T = 150 K, T = 300 K). Nevertheless, 
even at T = 150 K we can see non-vanishing Pt probability 
at the outside shell and some migration of Fe atoms towards 
the interior of the particle. This process seems fully consis¬ 
tent with our earlier observations drawn from the structures in 
Fig. [4] The high temperature in the T = 1000 K case helps to 
overcome the activation barriers and speeds up the reconstruc¬ 
tion. We expect that if we leave the lower temperature cases to 
evolve long enough they will reach similar shape of the distri¬ 
bution as the high temperature case. Unfortunately, this kind 
of time scale is prohibitively long to achieve - we estimate it 
will require at least one year worth of MD run for T = 150 K. 
Note also that after reconstruction some level of ordering is 
recovered in the high temperature system - the peaks are bet¬ 
ter defined and localised and actually narrower than for the 
low-temperature systems. In fact, one can recognise a clear 
gradient of decreasing noise and disorder of the distributions 


with increasing temperature in the right column of Fig. [5] In¬ 
terestingly, Fe atoms can be found close to the centre only at 
low temperatures, while they are systematically replaced by 
platinum atoms at higher temperatures. 

The results obtained for two types of particles can be easily 
understood if we take into account strong tendency of Pt atoms 
for occupying the outer layers of NPs. It prevents disordering 
of Pt-terminated particles, which keep their perfect multi-shell 
structure even at high temperatures. It agrees with the previ¬ 
ous theoretical studies 22 27 and the experimental observations 
of good thermal stability of icosahedral Fe-Pt particles even 
with much larger sizes. 20 In the thin films, the Pt-termination 
was deduced from the DFT calculations of the phonon spectra 
and the surface sensitive nuclear inelastic scattering measure¬ 
ment. 47 The simulated Fe-terminated surface showed strong 
deviation from the bulk behaviour. 

The crucial role in the stabilization of multi-twinned 
nanoalloys is played by the difference in atomic sizes of the 
two elements: larger size mismatch reduces the compression 
of the core. The MD simulations performed on the polyicosa- 
hedral core-shell Ag-Cu and Ag-Ni clusters revealed much 
higher melting temperatures than for pure Ag, Cu and Ni par¬ 
ticles. 4849 Even single impurities can stabilize the icosahe¬ 
dral Ag particles - the smaller the doping atom, the higher the 
melting temperature. 50 As shown in recent studies, the stabil¬ 
ity of Ni-Au NPs is also enhanced in the core-shell structure, 
which has much higher melting temperature than random al¬ 
loys. El 

The Fe 43 Pt 12 particle is very unstable due to the movement 
of Pt atoms towards the surface. This leads to significant 
changes of atomic positions and strong deformation of the 
icosahedral geometry. Apart from the segregation processes 
resulting from different sizes of constituent atoms, other amor- 
phisation mechanisms may also be of some importance here. 
The previous MD simulations revealed the amorphisation of 
mono-atomic Pt 55 icosahedral particles involving the rosette¬ 
like structural transformation, which forms a sixfold ring cen¬ 
tred around the fivefold vertices and breaks the 4 symme¬ 
try. 52 A preference for low-symmetry amorphous structures 
was found also in the Au nanoclusters. 53 

4 Disordered Fe-Pt nanoparticles 

The results described above concerned icosahedral multi-shell 
structures. These structures are inherently non-stoichiometric 
and thus quite different from the fragments of the bulk crys¬ 
tals. On the other hand, such highly ordered small particles are 
very unlikely to form naturally, although advanced experimen¬ 
tal techniques may be well capable of producing just such sys¬ 
tems. The synthesised icosahedral particles with larger sizes 
(^5 — 6 nm) and close to equiatomic compositions exhibit the 
shell-periodic structure with the Fe/Pt core and the Pt enriched 




outer shell. 20 For small systems such ordered core-shell struc¬ 
tures may be very unstable and we expect natural processes to 
form rather disordered particles with the composition of ele¬ 
ments determined by their relative concentrations in the envi¬ 
ronment and other factors such as chemical potentials, surface 
energies, bonding energies etc. 



Fig. 6 Small disordered particle at T = 300 K (a), larger disordered 
particle at T = 300 K (b) and at T = 500 K (c). All structures after 
30 ps thermalisation. Elements: Fe - dark/brown, Pt - light/grey. 


of the MD run) for both sizes of the particles are depicted in 
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As the simplest model we have decided to select two magic 
number (i.e. numbers from the partial sum: 1 + 12 + 42 + ...) 
particles with concentrations close to 50/50 (the magic num¬ 
bers are odd, thus the exact 50% concentration is unattain¬ 
able). The selected structures are Fe 6 Pt 7 and Fe 27 Pt 28 - both 
with slight iron deficiency - as we expect such particles to 
be more stable than their platinum-deficient counterparts (see 
above). We selected magic-number particles (13 and 55 atoms 
respectively) to make it possible for the systems to relax to 
closed-shell configurations which one can expect to be ener¬ 
getically favourable. 

The applied procedure was identical to the one used above 
in section [2^2] The only notable difference was a much longer 
thermalisation period required for the system to settle down 
close to the equilibrium - of the order of 10 ps. Due to a large 
size of the configuration space of disordered particle and a 
long time scale of slow relaxation processes, we do not expect 
the system to reach the true equilibrium state. Nevertheless, 
by examining the directions the system configuration drifts 
after thermalisation, we can determine the direction of these 
relaxation processes. This idea can be further cross-checked 
by application to very small particles. Due to much smaller 
configuration space and much lower barriers for collective 
movements, a small particle can reach its equilibrium in much 
shorter time. Running the MD simulation until full relaxation 
is actually feasible in this case. We ran this procedure for the 
smallest non-trivial, magic cluster: Fe 6 Pt 7 at T = 300 K and 
T = 500 K. The parameters of the simulation were the same 
as for larger systems, except that the total simulation length of 
30 ps was easy and inexpensive to achieve. Simultaneously, 
the small system stopped to drift in any significant way after 
only few picoseconds of simulation - indicating arrival at the 
equilibrium configuration. The resulting structures (at the end 


Fig. 7 Radial distribution of atomic positions from the mass centre 
(c.m.) of the system. The results for all temperatures are compared 
side-by-side for two investigated systems small (left) and large 
(right). The semitransparent colour of the Fe plot (orange) is used to 
show the Pt plot (blue) in the same range. 


We applied the same analysis techniques as described above 
for icosahedral particles. The resulting radial distribution 
functions are plotted in Fig. [7] The results for the small system 
indicate the same pattern we observed in the icosahedral parti¬ 
cles - the platinum atoms tend to move to the external layer of 
the particle and the iron atoms gather closer to the centre. The 
same pattern is less visible in the data for the larger particle 

- which is not as close to equilibrium as the small one - but 
the same trend is visible anyway. There is clear difference in 
the level of ordering exhibited by small versus larger particles, 
for the reasons explained above, but even for a larger system 
there are clear signs of the emerging structure in the form of 
obvious peaks in the RDF function and visible movement of 
platinum atoms to the outside of the particle. The details of the 
RDF distribution are subject to standard short-term statistical 
and thermal fluctuations, but the overall shape of the distribu¬ 
tions proved to be fairy robust and subject only to long time 
scale relaxation processes. We believe that if the system is 
left to evolve even longer, it will probably reach partially or¬ 
dered equilibrium with platinum atoms located mainly on the 
surface of the particle. Unfortunately, the expected time scale 

- measured at least in nanoseconds, perhaps even longer - is 
prohibitive for the pure DFT MD technique. This type of cal¬ 
culation would require a method which is at least few orders 
of magnitude more effective than pure DFT MD - e.g^neural 
network derived, multi-parameter effective potentials. 54 56 
























5 Conclusions 


In the presented work we have investigated the stability and 
dynamics of icosahedral nanoparticles constructed with two 
radically different (in terms of mass, atomic radius, bulk mod¬ 
ulus, magnetic moment) metals: iron and platinum. The op¬ 
timised icosahedral NPs with perfect layered structures have 
rather stiff platinum shells with inter-atomic NN distances 
well corresponding to the Pt bulk values and strongly stretched 
iron shells with NN distance elongated by 8 — 10% in compar¬ 
ison to the bulk. Even for iron-separated NPs, with significant 
dominance of Fe atoms, compression of Pt shells does not oc¬ 
cur and proves to be impossible. On the other hand, unstable 
iron shells drive the system to reconstruction and/or significant 
lowering of the melting temperature. 

The DFT molecular dynamics calculations carried out in a 
number of temperatures (T = 150 — 1000 K) revealed strong 
instability of the iron-terminated structure and significant ten¬ 
dency for platinum atoms to migrate into surface layer - even 
at low temperatures and platinum-deficient particles. The 
same tendency was observed in the very small (13 atoms) 
as well as larger disordered particles (55 atoms) with atomic 
concentrations close to the 50/50 ratio. Furthermore, the plat¬ 
inum capping layer seems to additionally stabilise the particle. 
In our calculations the platinum-terminated system proved to 
be stable up to T = 1000 K. On the other hand, the iron- 
terminated system shows substantial instability - leading to 
significant reconstruction of the structure and probably partial 
melting even at a low temperature (T = 150 K). 

The stabilizing effect of platinum atoms seem to be related 
to the geometric aspect of platinum enclosure of the NP. In 
the case of Fe 43 Pt 12 NP, platinum atoms have no significant 
influence on the stability of the NP because they only decorate 
the surface instead of forming complete enclosure. The case 
of 50/50 Fe-Pt NPs further indicates that it is rather geometry 
than concentration of platinum atoms that effects the stability 
of nanoparticles. 

The structural properties of small Fe-Pt NPs and their de¬ 
pendence on temperature discussed in the present work may 
be important for nano-technological applications. For in¬ 
stance, the catalytic activity of strongly distorted NPs could 
be enhanced comparing to pure Pt surfaces or larger Fe-Pt par¬ 
ticles. A larger distribution of Fe-Pt bond lengths as well as 
the presence of both types of atoms at the surface may lead 
to selective modification of electronic structure and catalytic 
activity of Pt atoms. It was observed that additional Co atoms 
decorating the surface of Pt nanocrystals lead to increased se¬ 
lectivity of hydrogenation reactions J 57 ^ 
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